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Abstract 

The influence of a random field induced by impurities, boundary irregularities etc. on the super- 
fluidity of a quasi-two-dimensional {2D) system of spatially indirect excitons in coupled quantum 
wells is studied. The interaction between excitons is taken into account in the ladder approxima- 
tion. The random field is allowed to be large compared to the dipole-dipole repulsion between 
excitons. The coherent potential approximation (CPA) allows us to derive the exciton Green's 
function for a wide range of the random field, and the CPA results are used in the weak-scattering 
limit, which results in the second-order Born approximation. The Green's function of the collec- 
tive excitations for the cases of (1) equal electron and hole masses and (2) the "heavy hole" limit 
are derived analytically. For quasi-two-dimensional excitonic systems, the density of the super- 
fluid component and the Kosterlitz-Thouless temperature of the superfluid phase transition are 
obtained, and are found to decrease as the random field increases. This puts constraints on the 
experimental efforts to observe excitonic superfluidity. 
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I. INTRODUCTION 



Superfluidity in a system of spatially indirect excitons in coupled quantum wells 
(CQW) has been predicted by Lozovik and Yudson/ and several subsequent theoretical 
studie o^i^i'^i^i^i'''i^i^i'^°i'^^i^^i^^i^^i^^ have suggested that this should be manifested as persistent 
electric currents, quasi- Josephson phenomena and unusual properties in strong magnetic 
fields. In the past ten years, a number of experimental studies have focused on observing 
these behaviors .^^i'^''i'^^i'^^i^°i^'^i^^i^^ The coupled quantum well system is conceptually simple: 
negative electrons are trapped in a two-dimensional plane, while an equal number of positive 
holes is trapped in a parallel plane a distance D away. One of the appeals of this system is 
that the electron and hole wavefunctions have very little overlap, so that the excitons can 
have very long lifetime (> 100 ns), and therefore they can be treated as metastable particles 
to which quasiequlibrium statistics apply. Also when D is large enough, the interactions 
between the excitons are entirely dipole-dipole repulsive. 

Much of the theory of excitonic systems has not taken into account the role of disorder, 
which is created by impurities and boundary irregularities of the quantum wells. In real 
experiments, however, disorder plays a very important role. Although the inhomogeneous 
broadening linewidth of typical GaAs-based samples has been improved from around 20 meV 
to less than 1 meV,— the disorder energy is still not much smaller than the exciton-exciton 
repulsion energy. At a typical exciton density of a few times 10^° cm~^, the interaction 
energy of the excitons is several meV. On the other hand, the typical disorder energy of 
1 meV is low compared to the typical exciton binding energy of 5 meV. Typical thermal 
energies at liquid helium temperatures are ksT = 0.2 — 2 meV. 

If the chemical potential of the interacting exciton system (controlled by the characteristic 
dipole-dipole interaction at fixed exciton density) is smaller than the characteristic disorder 
energy, the exciton system is expected to be localized in separate lakes in the minima of 
the random potential created by the disorder. The Bose condensation in this case is similar 
o that for Bose atoms in trap ^^i^^i^^i^'^'i^^i^^ (for the role of exciton tunneling, etc., see Ref. 
3ol|). In the opposite limit, one has a translationally invariant, extended exciton system 
with a random potential field. In the latter case, the critical properties and quasi-long non- 
diagonal order are similar to the case without disorder (see, e.g., Ref. jsil and references 
cited therein), but the disorder suppresses the Kosterlitz-Thouless transition temperature 
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and the superfluid density. 

Earlier studies of disorder in exciton systems included theory of the transport properties 
of direct and indirect excitons and magnetoexcitons in random fields,— the influence of var- 
ious random fields on excitonic and magnetoexcitonic absorption of light , and Anderson 
localization of excitons.-^ The effect of a weak random field on the collective properties and 

n 

superfluidity of excitons in nonuniform systems was analyzed in Ref . [36|| , including the cases 
of a dilute gas of three-dimensional (3D) excitons, two-dimensional (2D) excitons in a single 
quantum well, and indirect excitons in coupled quantum wells in a random field. In that 
work the random field was assumed to be much smaller than the exciton-exciton interaction. 
In two-dimensional systems, the excitonic interaction in the Bogoliubov approximation is 
valid only in non-physically low densities because of the divergence of the two-dimensional 
scattering amplitude in the Born approximation.™ Therefore, the ladder approximation 
must be used at low densities to treat properly the interaction between two-dimensional 
excitons.—"^ In Ref. , the contribution to the exciton Green's function from the interac- 
tion of the excitons with the random field was derived by perturbation theory, which limited 
the strength of the random field that could be studied. In this paper we study the case of a 
random field which is not necessarily small compared to the dipole-dipole repulsion between 
excitons. The coherent potential approximation (CPA) allows us to derive the 2D indirect 
exciton Green's function for a wide range of random field strengths, resulting in the second 
order Born approximation in the weak scattering limit (the second order Born approxima- 
tion Green's function for 3D direct excitons was obtained by Gevorkyan and Lozovib^^). 
We predict that in the low-temperature limit, the density Us of the superfluid component 
in CQW systems and the temperature of the superfluid transition (the Kosterlitz-Thouless 
temperature T^) are decreasing functions of the random field. 

A typical example of a two-dimensional system of weakly interacting bosons is a system 
of indirect excitons in coupled quantum wells (GaAs/AlGaAs) -^^'^'^'^^'^^'^'-''^^'^^ . Fluctuations 
of the thickness of a quantum well, which arise during the fabrication process, impurities 
in the system, and disorder in the alloy of the barriers can all lead to the appearance of 
a random field. Of these, spectral analysis of the exciton luminescence shows that alloy 
disorder, with a characteristic length scale short compared to the excitonic Bohr radius of 
around 100 A, plays the most important role. 

The paper is organized in the following way. In Sec. |n]the Green's function of the single 
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exciton in the random field is derived analytically in the coherent potential approximation 
(CPA) in the weak-scattering limit (resulting in the second order Born approximation) for 
the cases of (1) equal electron and hole masses and (2) the "heavy hole" limit. In Sec. 11111 
the Green's function and the the spectrum of collective excitations for a system of dilute 
indirect excitons in a random field are obtained within the ladder approximation. In Sec. II VI 
the dependencies of the density of the superfiuid component and the temperature of the 
Kosterlitz-Thouless transition on the random field are derived. In Sec. |3 we we present 
our conclusions and discuss possible experimental manifestations of superfluidity of indirect 
excitons in CQWs in a random field. 

II. THE GREEN'S FUNCTION OF THE SINGLE EXCITON IN THE RANDOM 
FIELD 

In our model, the random potential V{r) acting on an electron and hole is considered to 
be delta function correlated Gaussian noise, such that 



where d is the dimensionality of the space (for the spatially indirect exciton in CQWs d = 2). 
An electron is subjected to the potential Ve = aeV{r) and a hole to Vh = —ahV{r), where 
ae and at are constants. 

We consider the characteristic length of the random field potential I to be much lower 
than the average distance between excitons ~ 1/ ^/7m (/ ^ 1/ ^/^Tn, where n is the 
total exciton density). Therefore, in order to obtain the Green's function of the excitons 
with dipole-dipole repulsion in the random field, we first obtain the Green's function of a 
single exciton in the random field (not interacting with other excitons), and then apply 
perturbation theory with respect to the dipole-dipole repulsion between excitons, using the 
system of the non-interacting excitons as a reference system. 

The Green's function of the center of mass of the isolated exciton at T = in the 
momentum-frequency domain (G*^°^(p, cj)) in the random field in the coherent potential 
approximation (CPA) is given by^^ (here and below ^ = 1) 



{V{r)V{r'))=g5'ir-v') 



{V{r)) = 



(1) 



1 



(2) 



iu - eo{p) + /i + iQ{p,uj)' 
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where /i is the chemical potential of the system, and eo{p) = /2M is the spectrum of 
the center mass of the exciton in the "clean" system (M = me + nih is the mass of the 
exciton; mg and are the electron and hole masses, respectively). The function Q{^,uj) 
is determined by the effective random field acting on the center of mass of the exciton. 
For zero random field, Q{^^uj) — > 0. If gM ^ where Ei, is the binding energy of the 
indirect exciton (Eb ~ e^/tD at D > D being the distance between e and h wells; 

p{D) = {SaY^^D^^^ is the radius of the 2D indirect exciton at D ^ a, e is the dielectric 
constant, a = e/Arrie-he^ is the two-dimensional excitonic Bohr radius in the limit D ^ 0; 
e is the electron charge, and nie-h = 'merrih/ {jne + trih)), then the function Qijp^uj) in the 
coherent potential approximation is given by2^ 



g(p,o.) = ^/GW(q,c.)i?(|p-q|)^, 



(3) 



where 



S(p) = j d'^rB{R)e-'P^, (4) 

with R = Ri — R2; in the coordinate domain i?(Ri,R2) has the forn:^ (it will be shown 
below that -B(Ri, R2) depends only on R = Ri — R2) 
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where '^^{r) is the ground-state wave function of an exciton, corresponding to the relative 
motion of the electron and hole (r = rg — Vh)- For a two-dimensional indirect exciton with 
spatially separated electron and hole at large interwell distances D {D ^ a) the ground-state 
wave function is given by^ 
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exp 



7rp2(D) ^ \^ 2p2(D) J ■ 
Substituting foi^r) from Eq. © into Eq. (jSJ, and assuming = 2, we obtain 
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Using Eq. we obtain the Fourier transform of B{R) 

ByP) = ^ I ^ j ' = 

PTfn m.i. ^ m..- rvi. \ 



^(P) = T7^ ^^P — TT' ' rrih'^me] ah[ < Oe- (8) 

167r^ \ 16 y Vme/ 

Thus the CPA Green's function of the 2D indirect exciton is determined by the solution of 
the self-consistent equations Eqs. (j2I) and Q- In the weak-scattering limit ( q <^ /eDM ) 
we use the second-order Born approximation for Q similar to Refs. expanding Q 

(Eq. (jni) in a Taylor series to the first order in -B(|p — q|) (which is the first order in g), 
and we replace Eq. (jH)) by: 



(9) 



rUe = mh] 



Substituting B{p) from Eq. (jH)) into Eq. 0, we obtain for Q{p,Lj) 

nih^nie; ah[ < Oe, (10) 

where Jo{z) is a Bessel function of the first kind. The second-order Born Green's function of 
the single indirect exciton in the random field, G^^\p^ cu), is derived by substituting (5(p, oj) 
from Eq. jTUl) into Eq. (0). 



III. THE SPECTRUM OF COLLECTIVE EXCITATIONS 

At small densities n {np^ <^ 1), the system of indirect excitons at low temperatures is a 
two-dimensional weakly nonideal Bose gas with dipole moments d normal to wells {d ~ eD). 
The distinction between excitons and bosons manifests itself in exchange effects (see, e.g., 
Refs.Q.Q. andi^). These effects a.e suppressed fo. exe.tons witK spatially separated e 
and h m a, dilute system {np^ <^ 1) at large D {D ^ a), because at large D, the exchange 
interaction in the spatially separated system is suppressed relative to the e — h system in a 
single well due to the smallness of the tunneling exponent T ~ exp [-(D/2a)i/^] connected 
with the penetration of the dipole-dipole interaction through the barrier between the two 
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wellsi^. Hence, when D ^ a, exchange phenomena, connected with the distinction between 
excitons and bosons, can be neglected, and therefore, the system of indirect excitons in 
CQWs can be treated by diagram techniques employed for boson system^^. Two indirect 
excitons in a dilute system interact as U{R) = e^D"^ /{eB?), where R is the distance between 
exciton dipoles along quantum well planes. 

As mentioned above, in the two-dimensional case the contribution of the dipole inter- 
actions can be represented by the sum of the ladder diagrams given in Fig.l. The Bogoli- 
ubov approximation for a two-dimensional weakly interacting Bose gas is not valid due to 
the divergence of the two-dimensional scattering amplitude in the Born approximation^. 
Therefore, for np^ <^ 1, we take into account the direct dipole-dipole repulsion be- 
tween excitons within the framework of two-dimensional Bose gas theory in the ladder 
approximatioii^. Since the characteristic frequencies and momenta which give the great- 
est contribution to the single exciton Green's function G'*-°^(p, cu) in the ladder approxima- 



tion are^l ujeD/e^ ~ n/ \og[a^ /{nD"^)] < 1 and pp{D) ~ My'n/ log[aV(nD4))] < 1 (at 
nValP <C 1 and D ^ a), respectively, we approximate Q{p,uj) by Q{p = 0,a; = 0) (see 



Eq. m: 



Q(p,u;) = Q = ' = 

Q(p,u;) = Q = — — -, mh:^me] ah[ — < Og. (11) 



647r^ \meJ 

Note, that the replacing Qijp.uj) by constant Q(p = 0,u; = 0) does not constitute another 
additional approximation, but follows from Eq. (|l(Jj) at small frequencies and momenta, 
which is applicable for the dilute 2D dipole gas as mentioned above. We show below that 
the constant Q actually does not affect the ladder approximation vertex r(p, p', 0) for the 
weakly interacting Bose gas^. 

We pursue the 2D ladder approximation for the two-particle vertex r{p,p']P) (Fig.l) 
analogously to Ref. |37| 

Tip,p'; P) = f/(p - P ) + ^ / (0if^(P - + G^'^ - ^) (^2) 

where the arguments of the vertex and the Green's function are 3D momentum-frequency 
vectors (e.g., p = {p,uj}), and s is the level degeneracy (equal to 4 for excitons in GaAs 
quantum wells). Using the single exciton Green's function G^^\p, cu) determined by Eq. (j21) 
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with Q{p,u!) = Q from Eq. (fTT|l . we obtain the integral equation for F: 

A f/(p-q)r(q,p';P) 



r(p,p';P) = U{p-p') + s J 



{2tt 



|2 k£ 



M 



K 



/i 



2M 



nlo = nr(0,0;0), 



(13) 



where P = {P, fl}, and /i is the chemical potential of the system. 

The integral equation Eq. (jlHj) for the vertex can be solved analytically in the approxi- 
mation K, <^ ^/n.^'^ This inequality must be fulfilled simultaneously with the condition of low 



density np^i^D) ^ 1 (at ny/alP <^ 1 and D ^ a) which is necessary for the applicability 
of the ladder approximation. The solution of the integral equation for the vertex F of this 
system can be expressed through the solution of the equation for the scattering amplitude 
/o(k) = /(/t, k) of an isolated pair of interacting particles in a two-dimensional system with 
the repulsive potential U{R) = e^D"^ / {eR^): 

/o(-) = , yiLn^y (14) 



log 



The relation of the vertex r(p,p',0) (Fig.l.), which represents the effective interaction, to 
the two-dimensional scattering amplitude /o(p',p) is 
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El 

M 



F(K,p';0) 
(15) 



+ 2iQ 

where 5 — 0. In the first order in the scattering amplitude /o of the weakly interacting Bose- 
gas we neglect the second term in the r.h.s. of Eq. Therefore, the ladder approximation 
vertex F does not depend on the constant Q for the weakly interacting Bose gas to first order 
in the scattering amplitude. Hence, the vertex of the dipole repulsion for the rare exciton 
gas in the random field will be the same as for the clean system without random field. 

Here the characteristic momentum k, unlike in the three-dimensional system, is not equal 
to zero but rather is determined from the relation^i 

f2T^K\^/^ 



K = -Anfoin) 



(16) 



This is a specific feature of two-dimensional Bose system connected with the logarithmic 
divergence of the two-dimensional scattering amplitude at zero energy. A simple analytical 
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solution for the chemical potential can be obtained if nMe^D'^/e ^ 1. In this limit, the 
chemical potential ^ takes the form: 

2M 2Aflog(5;,;,J^)- ^ ' 

Since the only difference between the Green's function of the isolated exciton with and 
without the random field is the term iQ in Eq. all internal blocks of the ladder approx- 
imation diagrams for the system in the random field give the same self-energy as without 
the random field. This can be shown by repeating the procedure of the derivation of the 
2D ladder approximation self-energy for the "clean" system^J, by using integration measure 
d{uj + iQ) instead of duo. The difference in the Green's function will appear only in the 
external lines as the replacing of du by d{(jj + iQ). Introducing the Green's function of the 
Bose condensate and normal and anomalous Green's functions of the noncondensate analo- 
gously to Ref. 3^, we use the finite temperature Green's functions, replacing uj in the zero 
temperature Green functions by iuok, where ujk = 27rkT {k is an integer~; we set the Boltz- 
mann constant fc^ = 1). This procedure is valid because at small temperatures (T <^ /i) we 
assume the collective spectrum to be the same as the zero-temperature spectrunt^*^^^. For 
the 2D system the temperatures where superfluidity exists can be assumed small, because 
they are required to be below the Kosterlitz-Thouless temperature (we show below that 
Tc<fi). 

We therefore have the condensate Green's function D{p,iu!k) 

c'°'(P.^^>)- :'"!''"fn - (18) 

iujk + n + iQ 

where Uq is the density of Bose condensate. Since at small temperatures (n — nQ)/n ^ 
1, according to the ladder approximation's we use n below instead of hq. G{p,iujk) and 
F{p, iujk) are the normal and anomalous Green functions of the overcondensate: 

iujk + eo(p) + fJ' + iQ 



G{p,iuk) 



u;2 + e2(p)_2z(^-£ob))Q' 



where £o(p) is the spectrum of noninteracting excitons; the spectrum of interacting excitons 



F{p,iujk) = 2 , 2f \ — TT^T^' (1^) 

^k + ^ [P) - - eo{p))Q 



has the form e{p) = J (p'^/{2M) + ii^ — — (/x^ — Q^), and for small momenta p <ti fj, 



the excitation spectrum is acoustic e{p) = CsP, where Cs = y V f^'^ ~ Q'^/M is the velocity of 
sound. 
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From the expression for the spectrum e{p) we see that f or Q > fi the Bose condensate 
state becomes unstable, because the spectrum e{p) becomes imaginary. So, for Q > yU the 
random field destroys superfluidity, even at T = 0. This condition for the instability of 
the condensate itself is very approximate, because we have used the ladder approximation, 
which is valid only if almost all paricles are in the condensate {{n — nQ)/n 1), which is 
not the case for large random field contribution Q. Since at Q < fi the spectrum of the 
system is acoustic and satisfies the Landau criterium of superfiuidity, the system becomes 
a two-component liquid, consisting of the superfiuid and normal component in the presence 
of the random field even at T = 0^. As Q grows, the system undergoes a transition to 
the exciton glass state (for the a lattice model of bosons, this transition was considered in 



Ref. 



IV. THE KOSTERLITZ-THOULESS PHASE TRANSITION 

The density of the superfiuid component UsiT) can be obtained using the relation ns(T) = 
n — nn{T), where nn{T) is the density of the normal component. 

The density of the normal component nn{T), which is dissipated at the walls and impu- 
rities, can be calculated using the Kubo formula as the response of the total momentum to 
an external velocity^: 

'Im(n{iu;)y 



nn = - lini 



lUJ 



(20) 



where n(iu;) is the polarization operator with zero transferred momentum 



and J?^(p, iuJk) is the total single-particle Matsubara Green's function of an indirect exciton 

^(P> ^^'k) = D{p, io^fc) + G(p, iuj[). (22) 

The renormalization of the vertex by the interaction is neglected in the polarization operator 
Eq. ()2H1 . When the interaction is taken into account in the ladder approximation, a term 
which is small with respect to the parameter MT <^ 1 appears (F is the vertex in the ladder 
approximation). For a two-dimensional rarefied system of indirect excitons this parameter 
has the form An/ log [{e^ /{Sns'^nM^e^D^)] < 1. 
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We now substitute the Green's functions of the condensate Eq. (fTHj) and noncondensate 
Eq. (fT^ particles into Eq. (j221). Next, substituting the expression Eq. into Eqs. (PT|) 
and we have 



N 
M 



dp 2 ^o(p)^ 



(277)2- ^e^{p)^- ^^^^ 
Here is the total number of particles, and n° is the density of the normal component in 
a pure system with no impurities: 

1 



2M J {2tt 



dp 29no{p) 



de 



(24) 



oe{p)/T 



is the distribution of an ideal Bose gas of temperature excita- 



where no(p) = {e" 
tions. 

The first term in Eq. ()23|). which does not depend on Q, is the contribution to the normal 
component due to scattering of quasiparticles with an acoustic spectrum in an ordered 
system at T 7^ 0. In a two-dimensional system, 

27r ct(n,Q)M' 

where C,{z) is the Riemann zeta function (C(3) — 1.202). The second term in Eq. is the 
contribution to the normal component due to the interaction of the particles (excitons) with 
the random field. 



nl + s- 



nQ 



3C(3) T3 



+ s- 



nQ 



'2Mcl{n,Q) " 2n ciin,Q)M ' "2Mclin,Q)' ^^^^ 
The density of the superfiuid component is n<j = n — nn- From Eqs. and ^IF^ we can 
see that the random field decreases the density of the superfiuid component. 

In a 2D system, superfluidity appears below the Kosterlitz-Thouless transition temper- 
ature Tc = 7ms/{2M)^ where only coupled vortices are present. Using the expressions 
(|25p and (j26p for the density of the superfiuid component, we obtain an equation for the 
Kosterlitz-Thouless transition temperature Tc. Its solution is 



32/M^y 
\ 27 V vrn' J 



1/3 



^ 32 (MT^y 
^^27l,^J 



1/3- 



1 - 1 



c 

2V3' 



(27) 



Here T° is an auxiliary quantity, equal to the temperature at which the superfiuid density 
vanishes in the mean-field approximation (i.e., ns{T^) = 0), 



'27m'4M' 
. 3C(3) , 



1/3 



(28) 
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In Eqs. (I2Z|) and (I2H1), n' is 



(29) 



V. DISCUSSION 

The dependence of the Kosterhtz-Thouless transition temperature as a function of the 
total exciton density n for different Q, obtained from Eq. (j2Z|), is presented in Fig. 2. It can 
be seen in Fig. 2 that the random field decreases the Kosterlitz-Thouless transition temper- 
ature. This trend was pointed out before, for weak coupling with the random field, in Ref. 
3^. Figs. 3 and 4 show the dependence of the Kosterlitz-Thouless transition temperature 
on the random field parameter Q and the spatial separation between the electrons and holes. 

The results of the approximation used in Ref. 0| can be obtained from the normal 
density ?t,„, derived in the present work (Eq. (f^ ). as a first order in the expansion of 
in Eq. ()26|) respect to Q/fi, which corresponds to the case when the random field is weaker 
than the dipole-dipole repulsion and Q/fi<^ 1. (Operationally, this corresponds to replacing 



the speed of sound by its value at Q = 0, i.e., Cs yfj./M, in Eq. (j26|l .) For reahstic 
experimental parameters, the random field is not always smaller than the dipole-dipole 
repulsion. (E.g., in Ref. |42|, the luminescence linewdith due to inhomogeneous broadening 
at low density is approximately 2 meV; more recent GaAs structures'^ have inhomogeneously 
broadened linewidths closer to 1 meV.) Figs. 5 and 6 show that the approach used in the 
present work results in the Kosterlitz-Thouless temperature being smaller than the transition 
temperature obtained from the approximation used in Ref. , which is denoted here by PT 
(for perturbation theory; see below). The difference in the Kosterlitz-Thouless temperature 
between these two approaches increases when random field Q increases and exciton density 
n decreases (Figs. 5 and 6). The results of this comparison are reasonable, because the 



approximation used in Ref. j36| implies a first-order perturbation theory with respect to 
Q/yU for the Green's function, and the inequality Q/fi<^l holds for small Q and not very 
small densities n {fi increases as n increases). Note that in the present work the parameter 
Q / ji\s not required to be small. 

This work shows that although the random field depletes the condensate, Kosterlitz- 
Thouless superfiuidity is still possible in a system of spatially indirect excitons. However, 
for realistic experimental parameters, the possibility of superfiuidity is marginal. For the 
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structure parameters used in the plots of Figs. 2-4, with effective D = 15 nm, the binding 
energy of the indirect excitons at high electric field is 3.5 meV and the excitonic Bohr 
radius is around 150 A. This implies that na"^ = 1, the point at which the exciton gas 
becomes plasma-like, at a density of around 5 x 10^^ cm^^. If the exciton density is to 
be well below this value, say an order of magnitude, then as seen in Fig. 2, the disorder 
parameter Q must be around 0.1-0.2 meV. By comparison, the half-width, half maximum 
of the exciton luminescence line at low temperature and low density in recent samples, 
due to inhomogeneous broadening, is around 0.5 meV— . It is not clear exactly how the 
inhomogeneous line width and the disorder factor Q are related. The line width is likely 
sensitive to disorder with length scales short compared to the excitonic Bohr radius, while the 
Q parameter used here is assumed to be a measure of disorder on length scales long compared 
to the excitonic Bohr radius. Nevertheless, there is still a problem with Kosterlitz-Thouless 
superfluidity at low density, no matter what the temperature. 

The story may change when the excitons are trapped in a confining potential, in which 
case true Bose condensation is possible in the Stringari-Pitaevskii limit. However, these re- 
sults indicate that the role of disorder must be taken into account for any realistic treatment 
of excitonic condensation. 
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Captures to Figures (1-6) 



Fig.l. Diagrammatic representation of the equation for the vertex T in tiie momenta- 
frequency representation (P,f2). 

Fig.2. Dependence of temperature of Kosterhtz-Thouless transition = Tc{n) (in units 
of K; for GaAs/AlGaAs: M = 0.24mo, where mo is the free electron mass; e = 13) on 
the exciton density n (in units of cm~^) at the interwell distance D = 15 nm, for different 
random fields Q (in units of meV): Q = - sohd curve (straight hue); Q = 0.1 meV - 
dotted curve; Q = 0.2 meV - dashed curve; Q = 0.3 meV - long-dashed curve; Q = 0.4 meV 
- dashed-dotted curve; Q = 0.5 meV - sohd curve (bottom right hand corner). 

Fig.3. Dependence of temperature of Kosterhtz-Thouless transition = Tc{Q) (in units 
of K; for GaAs/AlGaAs: M = 0.24mo; e = 13) on the random field Q (in units of meV) 
at the interwell distance D = 15 nm, for different exciton densities n: n = 5 x 10^^ cm~^ - 
solid curve; n = 1.0 x 10^^ cm~^ - dotted curve; n = 3.0 x 10^^ cm~^ - dashed curve. 

Fig.4. Dependence of temperature of Kosterhtz-Thouless transition Tc = Tc{D) (in units 
of K; for GaAs/AlGaAs: M = 0.24mo; e = 13) on the interwell distance D (in units of nm) 
at the exciton density n = 1 x 10^^ cm'"^ , for different random fields Q: (5 = 0- solid curve; 
Q = 0.3 meV - dotted curve; Q = 0.5 meV - dashed curve; Q = 0.6 meV - dashed-dotted 
curve. 

Fig.5. Temperature dependence of Kosterhtz-Thouless transition = Tc{n) based on 
CPA and the perturbation theory (PT) approximation used in Ref. 120] (in units of K; for 
GaAs/AlGaAs: M = 0.24mo; e = 13; mg is a mass of electron) on the exciton density n (in 
units cm~^) at the interwell distance D = 15 nm, for different random fields Q (in units of 
meV): Q = 0.2 meV - sohd curve, full CPA vs. dotted curve, PT; Q = 0.3 meV - dashed 
curve, full CPA vs. long-dashed curve, PT. 

Fig.6. Temperature dependence of Kosterhtz-Thouless transition Tr = Tc{Q) based on 
CPA and the perturbation theory (PT) approximation used in Ref. j36| (in units of K; 
for GaAs/AlGaAs: M = 0.24mo; e = 13) on the random field Q (in units of meV) at the 
interwell distance D = 15 nm, for different exciton densities n: n = 5 x 10^° cm~^ - solid 
curve, full CPA vs. dotted curve, PT; n = 1 x 10^^ cm~^ - dashed curve, full CPA vs. 
long-dashed curve, PT. 
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Figure 6 
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